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ABSTRACT 

We report on VLBI observations of the Type I b/c supernova 2001em, three years after the ex- 
plosion. It had been suggested that SN 2001em might be a jet-driven gamma ray burst (GRB), 
with the jet oriented far from the line of sight so that the GRB would not be visible from earth. 
To test this conjecture, we determined the size of SN 2001em. It is only marginally resolved at our 
resolution of ~0.9 mas. The 3er upper limit on the major axis angular size of the radio source was 
0.59 mas (FWHM of an elliptical Gaussian), corresponding to a one-sided apparent expansion velocity 
of 70,000 km s _1 for a distance of 80 Mpc. No low-brightness jet was seen in our image to a level 
of 4% of the peak brightness. If we assume instead a spherical shell geometry typical of a supernova, 
we find the angular radius of SN 2001em was 0.17^q'°q mas, implying an isotropic expansion velocity 

of 20, 000^2°ooo k m s_1 ' which is comparable to the expansion velocities of supernova shells. Our 
observations, therefore, are not consistent with a relativistically expanding radio source in SN 2001em, 
but are instead consistent with a supernova shell origin of SN 2001 's radio emission. 
Subject headings: (stars:) supernovae: individual (SN2001em), radio continuum: general, gamma rays: 
bursts 
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1. INTRODUCTION 



S upernova 2001em was d iscovered on 2001 September 
15 llPapenko va et al l 20011). at a distance of ~ 80 Mpc 
(e.g.. lSoderberg et al Il2004bt iHuchra et alJll992L assum- 



ing H = 70 km s" 1 Mpc" 1 ), in the galaxy UGC 11794 
(NGC 7112). It was classified as Type Ib/c on the ba- 
sis of an early optical spectrum ijFilipnenko fc Chornockl 
I2001D . A recent spectrum, however, was unusual for a 
Type Ib/c supernova, being dominated by a broad Ha 
line (FWHM of 1800 km s _1 ) suggestiv e of strong inter- 
action with the circumstellar medium l|Soderberg et all 
l2004b|). SN 2001em was r ecently also detected in the ra- 
dio ijStockdale et al.ll2004[K with an 8.4 GHz flux density 
of 1.15 mJy on 2003 October 17. It is unusually radio 
luminous, more so than any other observed Type Ib/c su- 
pernova at its age. Furthermore, SN 2001em's flux den- 
sity rose to 1.48 mJy by 2004 January 31, that is ~ 1000 d 
after the explosion, which is also unusual behavior for a 
Type Ib/c SN, whose radio light curves usually peak less 
than 1 year after the explosion. In addition, SN 2001em 
had a radio spectral index of a ~ —0.37 (S v oc v a ) be- 
tween 5 and 15 GHz, which is flatter than is usually seen 
in radio supernovae older than one year. 

The unusual ch aracteristics of SN 2001em prom pted 
the suggestion by Granot & R amirez-Ruizl l)2004l) that 
it might be associated with a gamma ray burst (GRB). 
There is mounting evidence that at least some GRBs are 
associated with supernovae. In particular, two GRBs 
have now been reliably associated with o bserved super- 
novae : GRB 980425 with SN 1998bw llGalama et al l 
1998) and GRB 030329 with SN 2003dh ifstanek et all 
2003; iHiorth et al.l 1200*3^ . both of which were of Type 
Ib/c. 

Most current GRB models rely on the collapse of 
a stellar mass object to a black hole, which powers 
relativistic jets (see e.g ., iMacFadven fc Wooselvl Il99fl 
Zhan g fe Meszarosll2004]) . GRB jets are thought to ini- 
tially be highly relativistic and have consequently highly 



beamed emission. GRB afterglows, particularly in the 
radio, are seen as the jets decelerate and their emission 
becomes more isotropic. If the jet model for GRBs is 
correct, it follows that only a fraction of the jets are ori- 
ented near the line of sight and are visible from earth 
as GRBs. The remainder will not produce a GRB event 
visible from earth. The more isotropic afterglow emis- 
sion, however, in particular radio emission from the jet, 
may be detectable. One might therefore hope to detect 
a signature of a GRB i n some Type I b/c supernovae. 
iSoderberg et alJ l)2004a|) have estimated the fraction of 
GRB events visible from earth as b eing < 6% (see also 
iBerger et af1l2003t iFrail et aUkOOlfl . Radio searches for 
GRB signatures in Type Ib/c sup ernovae were carrie d 
out bv lSoderberg et alJ l|2004aD and lBerger et alJ l|2003ft . 
but none have been detected to date. It was also pre- 
dicted that some GRB jets, especially ones not aligned 
near the line of sight, might be resol vable with very- 
long-baseline interfer ometry (VLBI) l|Paczvnskil 120011: 
IGranot fc Loebll2003l) . Recently, Taylor et al. (2004; see 
also Taylor et al. 2005) resolved GRB 030329 using VLBI 
observations, and showed that the radio source was ex- 
panding at an average apparent velocity of ~ 4 c, consis- 
t ent with the expectations . 

IGranot fc Ramirez-Ruizl |2004) showed that the char- 
acteristics of SN 2001em's radio emission, in particu- 
lar the delayed rise in flux density coupled with a ra- 
dio spectral index of a ~ —0.4, are difficult to recon- 
cile with the circumstellar interaction model of super- 
nova radio emission, but are consistent with the ex- 
pectations for an off-axis relativistic jet. Such a jet 
would be expected to become non-relativistic after a few 
years, with the radio emission re aching its maximum 
near the time of thi s transition l)Granot fc Saril 120021 
IGranot fc Loebll2003(K This behavior is consistent with 
the radio lightcurve of SN 200 lem. The recent rise in 
X-ray flux from SN 2001em l|Poolev fc Lewinl 12004ft is 
also consistent with the possibility that SN 200 lem was 
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an off-axis GRB. Could SN 2001em harbor a GRB jet 
transverse to the line of sight? 

At 80 Mpc, the length of a transverse jet which 
had expanded relativistically since the explosion date of 
SN 2001em (~ 2001 September 15) would correspond 
to an angular size of ~ 2 mas. The angular size of 
an ordinary supernova, by contrast, is expected to be 
< 0.2 mas, since supernova s hells expand with velocities 
of < 20, 000 km s" 1 (see e.g.. lBietenhoTzll2005T) . 

A clear observational distinction between radio emis- 
sion from a relativistic GRB jet and that from a super- 
nova shell is thus possible. In this letter, we describe 
VLBI observations of SN 2001em that we carried out in 
order to determine its size, and address the question of 
whether the radio emission is due to a relativistic GRB 
jet or a supernova. A n earlier attempt was made by 
iStockdale et alJ l)2005|) to determine a size with VLBI, 
who obtained an upper limit of 1.9 x 0.8 mas, correspond- 
ing to a non-conclusive upper limit of ~ 260, 000 km s _1 
or 0.9 c at 80 Mpc on one-sided expansion. 1 

2. OBSERVATIONS 

Because of the low flux density of only ~ 1.5 mJy, 
we used the "High Sensitivity Array" , which consists of 
NRAO's 2 Very Long Baseline Array (VLBA; each of 25 m 
diameter), phased Very Large Array, (VLA; 130 m equiv- 
alent diameter), and Robert C. Byrd telescopes (GBT; 
~105 m diameter), in addition to the Arecibo Radio Tele- 
scope 3 (Ar; 305 m diameter) and the Effelsberg Radio 
Telescope 4 (Eb; 100 m diameter). We phase-referenced 
our VLBI observations to obtain both reliable starting 
models for possible self-calibration and accurate astro- 
metric positions. Our primary phase-reference source 
was J2145+11, which is 1.4° away on the sky. How- 
ever, we spent two periods of ~20 min phase-referencing 
to J2139+14, (2.1° away) which is a "defining source" 
of the Internat ional Celestial Reference Frame (ICRF; 
iMa et akl 119981) in order to determine coordinates for 
SN 2001em tied to the ICRF. 

The VLBI observations were carried out on 2004 
November 22, with a total time of 12 hrs. We used a fre- 
quency of 8.4 GHz with a bandwidth of 32 MHz, and we 
recorded both senses of circular polarization. As usual, a 
hydrogen maser was used as a time and frequency stan- 
dard at each telescope. The data were recorded with 
the VLBA or the MKIV VLBI systems with a sampling 
rate of 256 Mbits s _1 , and correlated with the NRAO 
VLBA processor at Socorro. The analysis was carried 
out with NRAO's Astronomical Image Processing Sys- 
tem (AIPS). The initial flux density calibration was 
done through measurements of the system temperature 

1 IStockdale etall 120051) cite a limit on two-sided expansion of 
150,000 km s — 1 . It is likely, however, that the jet emission domi- 
nates the counter-jet emission, and hence we cite the more conser- 
vative velocity limit on one-sided expansion. We further convert 
their limit, which was derived using a distance of 90 Mpc, to our 
assumed distance of 80 Mpc. 

2 The National Radio Astronomy Observatory, NRAO, is oper- 
ated under license by Associated Universities, Inc., under cooper- 
ative agreement with National Science Foundation, NSF. 

3 The Arecibo Observatory is part of the National Astronomy 
and Ionosphere Center, which is operated by Cornell University 
under a cooperative agreement with the NSF. 

4 The 100 m telescope at Effelsberg is operated by the Max- 
Planck-Institut fur Radioastronomie in Bonn, Germany. 



(T sys at each telescope, and then improved through self- 
calibration of the phase- reference sources J2 145+11 and 
J2139+14. We used a cycle time of ~250 s, with -170 s 
spent on SN 2001em. The data from the Kitt Peak 
VLBA station were of very poor quality and not used 
in the final analysis. Finally, in addition to using the 
phased VLA, which was in the A configuration, as an el- 
ement in our VLBI array, we also used it as a stand-alone 
interferometer to determine the integrated flux density of 
SN 2001em. 

3. RESULTS 

Using the VLA as a stand-alone interferometer, we 
found the integrated flux density of SN 2001em on 2004 
November 22 to be 1.50 ± 0.10 mJy at 8.4 GHz, consis- 
tent with the ear lier value of 1.48 ±0.05 mJy obtained on 
2004 January bv IStockdale et aD l)2004[) . and suggesting 
that the flux density at 8.4 GHz is still increasing or at 
least staying relatively constant. 

In Figure 1 we show the VLBI image of SN 200 lem. 
No structure is apparent, and the source is at best 
marginally resolved at our resolution of 1.3 x 0.6 mas 
at p. a. —4° (FWHM of a Gaussian; p. a. is east of 
north). In particular, no s tructure large r than 1.3 mas 
or 1.6xl0 18 cm is visible. iTavlor et all l)2004|) found a 
radio component in GRB 030329 which, if it is related 
to GRB 030329, suggests an apparent expansion veloc- 
ity of —19 c. We made large images of 0.2" x 0.2", to 
search for such components, but found none: the bright- 
est points except for SN 200 lem were no brighter than 
expected from the noise, or < 0.1 mJy/bm. In addition, 
our phase-referenced VLBI observations give an accurate 
position for SN 2001em of R.A. = 21 h 42 m 23*60936(2), 
deck = 12° 29' 507 3003(3) (J2000; the figure in brackets 
gives the estimated uncertainty in the last digit). 5 

To get a more accurate determination of the size, we 
turned to model fitting. This process involves the least- 
squares fitting of geometrical models to the VLBI vis- 
ibility data in the Fourier transfor m plane, and is de- 
scri bed in more detail in, e.g., [Bietenholz et al. ( 2001) 
and lBartel et alJ ((2002) . It has the advantage of giving 
results independent of the convolving beam. 

Since the source is at best marginally resolved, we must 
choose the model a priori. We use two different models, 
the first being an elliptical Gaussian, which we choose to 
model a possibly elongated source such as a jet, and the 
second being a spherical shell model appropriate for a 
supernova. In particular, our shell model consists of the 
projection of a spherically symmetrical, optically thin 
shell of emission, with the shell width being 20% of the 
outer radius. Such a shell model is suggested by obser- 
vations of SN 1993.1 llBietenholz et al.ll2001tlBartel et alJ 
I2002t iBietenholz et all 1200.1 see also Marcaide et al. 

5 This position is that of the radio brightness peak. It 
was determined with respect to the phase calibrator source 
J2139+14, for which we assumed the ICRF position of R.A. = 
21 h 3 9 m If 30927(2), deck = 14° 23' 35'.' 9920(3) C .12000: IMa et alJ 
1998). Our estimated uncertainties are those from our phase- 
referencing analysis (150 /^as), added in quadrature to the listed 
uncertainties of the ICRF position for J2139+14 (~300 /xas). 
We note also that the differential position to our other reference 
source, J2145+11, is somewhat more accurate due to the higher 
signal-to-noise ratio and the somewhat smaller angular distance to 
SN 2001em. However, since the uncertainty of J2145+ll's position 
is 10 mas, the absolute position derived from it is less accurate. 
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Fig. 1.— A VLBI image of SN 2001em on 2004 November 22. 
The peak brightness was 1.46 mjy bm -1 and the rms background 
was 17 fijy bm" 1 . The contours are drawn at —4, 4, 10, 20, 30, 
50, 70 and 90% of the peak brightness, with the lowest contour 
being at 3<r, and the 50% contour being emphasized. The FWHM 
size of the Gaussian restoring beam, which can be compared to the 
50% contour of the source, was 1.30 X 0.62 mas at p. a. -4°, and 
is indicated at lower left. The greyscale is labeled in mjy bm -1 . 
The origin is at R.A. = 21 h 42 m 23 s 60936, deck = 12° 29' 50'.'3003 
(J2000). North is up and east to the left. 



1995, who find a somewhat thicker shell). 

The statistical uncertainties can be determined in the 
fitting process, since those of the visibility measurements 
are known. The visibility data, however, may suffer from 
additional uncertainties due to possible calibration er- 
rors, which are not statistically independent. We thus 
use the square root of the visibility weights in our fits, 
which results in a more robust fit less dominated by a 
few very sensitive baselines. The fit size is probably 
most sensitive to the amplitude gain calibration, since 
for a marginally resolved source like SN 2001em, the fit 
source size is correlated with the individual antenna am- 
plitude gains. Our antenna amplitude gains were de- 
rived first from the T sys measurements at each telescope 
and then refined by self-calibration using the reference 
sources. They should be un-biased, since they were de- 
rived without reference to SN 2001em itself, and accurate 
to < 5%. Our conservative uncertainties are obtained by 
letting the antenna amplitude gains be free parameters 
in the fit, which allows for possible correlations between 
them and the fit size. Our best-fit size, however, is the 
unbiased one arrived at by fixing the antenna gains at the 
values derived from the reference sources. This best-fit 
size is consistent with that derived by freeing the antenna 
gains. 

Using an elliptical Gaussian model appropriate for a 
jet, we obtained a FWHM major axis of 0.23lg mas; 
the 3cr upper limit on the FWHM major axis size is thus 
0.59 mas. The nominal p. a. of the major axis was ~ 97°; 



however we find that neither the ellipticity nor the p. a. 
is reliably determined by our data. 

For the spherical shell model appropriate for a super- 
nova, we find an outer angular radius of 0.17^g'5o mas, 
with a corresponding 3er upper limit of 0.35 mas. We 
note that other models with a geometry different from a 
spherical shell will give different outer radii, but for likely 
geometries, the differences are < 25%. For example, a 
uniform sphere gives a radius ~ 16% larger, whereas a 
very thin shell gives a radius ~ 10% smaller (see e.g., 
iBartel et al1l2002UMarscherlll985|) . 

4. DISCUSSION 

We have made a VLBI image of SN 2001em, a su- 
pernova with an unusually delayed onset of radio emis- 
sion, which was thought possibly to be due to a rela- 
tivistic jet not aligned with the line of sight. Such a jet 
would link SN 2001em to a GRB. The radio morphology 
of a relativistic j et wo uld be expected to be elongated. 
lOranot fc Loebl l)2003l) showed that the radio emission 
from such a jet would peak near the time when the jet 
expansion becomes non-relativistic. Thus up to the time 
of the radio peak, a transverse jet should be seen to ex- 
pand at an apparent speed of ~ c, with the apparent 
expansion speed being superluminal for a jet directed to- 
ward us, and sub- luminal for a jet directed away. Since 
the radio flux density has not yet started to decline sig- 
nificantly, we would expect such a jet in SN 2001em to 
still be in its relativistic expansion phase, and thus the 
brighter jet directed toward us should expand with an 
apparent speed of at least c. In the case that the jet in- 
clination angle is near 90° , we would expect the jet and 
counterjet to be of nearly equal radio brightness, and 
an end-to-end apparent expansion velocity of ~ 2 c. At 
80 Mpc, this implies a current angular size of > 2 mas 
for the approaching jet. 

At epoch 2004.9, we found a 3er upper limit to the 
FWHM major axis size of SN 2001em of 0.59 mas, corre- 
sponding to a size of < 7.1 x 10 17 cm. Assuming one-sided 
expansion, this corresponds to a projected average veloc- 
ity since the explosion of < 70,000 km s^ 1 or < 0.23 c. 
In the case of a two-sided jet, the velocity of each jet 
would be < 35, 000 km s -1 . Since no extended emission 
is visible in our image, any source expanding with a pro- 
jected velocity of > 0.5 c can be responsible for at most 
10% of the radio flux density. 

This upper limit on the expansion velocity implies that 
the bulk of the radio emission cannot be due to a rela- 
tivistic jet associated with SN 2001em. We note that a 
relativistic jet aligned near the line of sight is unlikely, as 
it would be seen to expand super-luminally. If there is a 
radio-bright jet, it must have ceased relativistic expan- 
sion at an age of < 0.5 yr. We cannot of course exclude 
the possibility of the body of the jet being faint in the ra- 
dio, with either the core or a terminal hot-spot being re- 
sponsible for the bulk of the radio emission. We note that 
future astrometric measurements, using the same refer- 
ence sources we used above, could allow a determination 
of the proper motion. In the case of a bright terminal 
hot-spot moving at > 40, 000 km s _1 , the proper motion 
would be > 0.1 mas yr -1 , which is easily detectable in a 
few years. 

The upper limit on the expansion velocity of 
SN 2001em is, on the other hand, consistent with those 
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usually seen in radio supernovae, in which the radio 
emission originates from the interaction between the su- 
pernova shell, which expands with velocities of order 
f 0,000 km s _1 , and the circumstellar medium. If we 
assume a spherical shell geometry for SN 2001em, we 
find its outer angular radius to be 0.17^q'^q mas, cor- 
responding to a radius of 2.0l°'2xl0 17 cm, and an ex- 
pansion velocity of 20, OOOl^ooo ^ m s_1 ( as & bove this 
value represents the average expansion velocity since the 
supernova explosion). This expansion velocity is con- 
sistent with tho se seen in other radio supernovae (e.g., 
lBietenholzll2005lL 

If the radio emission of SN 2001em is due to an ex- 
panding supernova shell, then we are left with the puzzle 
of its delayed onset, and also of the spectral index which 
is somewhat flatter than usually observed in supernova 
shell emission. The spectral index is, however, consistent 
with what is observed for pulsar nebulae. It is perhaps 



possible then, that the radio emission in SN 2001em is 
due to the turn-on of a pulsar nebula. As such it would 
be ~ 5000 times as radio-luminous as the 1000-year old 
Crab Nebula is at present. Such a high luminosity is not 
excluded for a y oung pulsar nebula ijBandiera et all!984: 
IChevalieri2 005) . We note that a new component with an 
inverted radio spectrum was rece ntly seen in SN 1986J , 
almo st 20 yr after the explosion ijBietenholz et al.ll2004l 
2002). This component, which is 200 times as radio- 
luminous as the Crab Nebula, may be a young pulsar 
nebula. Future monitoring of SN 200 lem will be neces- 
sary to determine whether the radio emission is due to a 
supernova shell, a pulsar nebula, or something else. 



We thank M. Rupen for useful comments on the 
manuscript. Research at York University was partly sup- 
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